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ABSTRACT 

Of the more than 150000 targets followed by the Kepler Mission, about 10% were selected as red giants. 
Due to their high scientific value, in particular for Galaxy population studies and stellar structure and evolu- 
tion, their Kepler light curves were made public in late 20 11. More than 13 000 (over 85%) of these stars show 
intrinsic flux variability caused by solar-like oscillations making them ideal for large scale asteroseismic inves- 
tigations. We automatically extracted individual frequencies and measured the period spacings of the dipole 
modes in nearly every red giant. These measurements naturally classify the stars into various populations, such 
as the red giant branch, the low-mass (M/Mq < 1.8) helium-core-burning red clump, and the higher-mass 
(M/Mq > 1.8) secondary clump. The period spacings also reveal that a large fraction of the stars show ro- 
tationally induced frequency splittings. This sample of stars will undoubtedly provide an extremely valuable 
source for studying the stellar population in the direction of the Kepler field, in particular when combined with 
complementary spectroscopic surveys. 

Subject headings: stars: interiors — stars: oscillations — techniques: photometric 



1. nvlTRODUCTION 

We have seen dramatic progress in the asteroseismology of 
red giant stars in recent years, driven by the quantum leap 
in data qualit y and quantity arising from t he space missions 
CoRoT (Baglinetal. 2006) and Kepler fcoch et alj 120101: 
iGilUland et al. 2010) . We can now start using large samples of 
stars to study stellar evolution in different populations within 
the Galaxy to a degree that has not been possible before. 

For about a decade it has been evident that red giants 
show oscillatio ns excited by n ear-surface convection like in 
the Sun (Fr andsen et al.ll2002h . known as solar-like oscilla- 
tions. This sparked hope that one could obtain rich infor- 
mation about the interior structure of red giant stars, analo- 
gous t o what global helio seismology had taught us about the 
Sun ( Christensen-Da lsgaard 2002). However, only recently 
was it demonstrated that the frequency spect ra of red giants 
comprised both radial and non radial modes (Ide Ridder et al.l 
I2OO91) . whose frequencies, f„ largely followed the asymp- 
totic rel ation for high-order, n, and low-degree, I, acou stic (p) 
modes (IVandakurovl 1 9671: iTassoulll 1 9801: iGoughll 1 9861) as for 
the Sun, 

v„^ic:i Au{n + l/2 + e)-Siyo,i. (1) 

Here, e is the offset from zero of the fundamental radial mode 
in units of the large separation, Ai/, which is the frequency 
shift of consecutive overtone modes of the same degree, and 
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Svqj is the small separation of non-rad ial modes rel ative to 
radial modes following th e notation o f (Bedding 2 01 1 | ). It was 
subsequently shown (Bedding et al.li2010: Huber et al.l 120101 : 
iMosser et al.l201 lai that the frequency spectra of the p modes 
in red giants are essentially similar under a simple scaling, 
which implies that the stars are largely homologous. 

In addition to the regular, and largely scalable, frequency 
pattern of the p modes, red giants show a more complex pat- 
tern arising from large numbers of mix ed modes, most pro - 
nounced for dipole, I = 1, modes ( Bedding et al.l l2010h . 
These mixed modes are the result of coupling between the 
many dipole gravity (g) modes in the core, which are approx- 
imately equally (and closely) spaced in period, and the more 
widely spaced dipole p modes in the env elope, which are ap - 
proximately equally spaced in frequency (iDupret et al.ll2009T) . 
While the coupling between p- and g-modes shifts the fre- 
quencies from their equal spacing, a process known as mode 
bumping, the period spacing of the mixed modes still is a 
prox y for the intrinsic period spacing of the g modes in the 
core dBeck et al.ll201 1 ). A detailed study of the difference be- 
tween the period spacing of uncou pled g-modes and th at of 
the mixed modes was conducted by iMosser et al.l (l2012ah . 

Due to their partial g-mode nature, the mixe d modes offer a 
'windo w' into the structure of the stellar cores. iBedding et al.l 
(I2OI ll) showed that their period spacing distinguish red gi- 
ant branch stars, burning only hydrogen in the shell, from red 
clump stars that also bu rn helium in their cores. More re- 
cently, iBeck et al] (1201 2h demonstrated that the mixed modes 
probe the internal differential rotation of red giants, enabling 
investigations into the transp ort of angular momen tum as stars 
'climb' the red giant branch (IMosser et al.ll2012bl) . 

The potential for obtaining new results on red giants was 
recently boosted by the public release of over two years of 
Kepler data for about 15 000 stars classified as giants in the 
Kepler Input Catalog - providing an order of magnitude more 
stars with such long time series than has previously been 
available through the Kepler Asteroseismic Science Consor- 
tium. Initial analyses were performed by Hekker et al. (201 ll) 
who measured Aiy and the location of the maximum oscilla- 
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tion power, J^max, of roughly 1 1 000 of these red giants based 
on the first 43 days of Kepler data. 

Here, we report a search for oscillations in this much longer 
data set in order to perform 'ensemble peak-bagging' - ex- 
tracting oscillation frequencies on large numbers of stars. 
From these we measure characteristic separations of the 
modes in frequency and period, which we analyze to iden- 
tify different stellar populations and demonstrate the potential 
for follow-up investigations of this extremely valuable stellar 
sample. 

2. DATA ANALYSIS 

Our initial data set comprised photometric time series of 
15 261 stars observed by Kepler from 2009 May 2 to 2011 
March 14 (observing quarters 0-8) in the long-cadence mode 
{5t = 29.4244 minutes), with a median time span of just be- 
low tw o years. Most of these stars are expected to be red 
giants dHekker et al.l l201 1). The large number of stars in our 
sample calls for an automated approach to the analysis, which 
we describe below. 

2.1. Measuring global seismic properties 

We have used simple aperture photometry (SAP) data 
for our analysis (iJenkins et al.l 12010 ), downloaded from the 
MAST database'. Jumps in each time series (e.g. across 
quarterly gaps) were corrected by fitting and removing a lin- 
ear regression to 5-day light curve segments before and after 
each gap. Remaining slow instrumental trends were then re- 
moved by subtracting a smoothed version of the light curve 
that was calculated by applying a boxcar filter with a width 
of 10 days. The high-pass-filtered light curves were then used 
to calculate a freque ncy spectrum. Thi s was analyzed with an 
automated pipeline (iHuber et al.ll2009l ) to detect the presence 
of solar-like oscillations. The results from the pipeline were 
verified by visual inspection of its graphical output of each 
star, and a few percent of the stars were discarded in this pro- 
cess because the pipeline results were deemed incorrect. Os- 
cillations were detected in 13 412 red giant stars. We found 
a similar number (131 82) when u sing the method described 
in lMosser & Appourchauxl(l2009l) . Hence, for ^ 2 000 stars, 
predominantly low and high log g giants, we 'missed' the os- 
cillations because of the high-pass filter and the limits set by 
the Nyquist frequency of the data, respectively. For each star, 
the analysis provided a frequency spectrum corrected for the 
stellar granulation signal, as well as measurement of I'max and 
Aj/ (FigurelUa)). 

2.2. Extracting oscillation mode frequencies 

Using an approach similar to iBedding et al.l (1201 1 *). we 
smoothed the frequency spectra and measured the frequency 
of each peak reaching above three times the noise level within 
±5 Ai/ of j/jnax- A smoothing width of about 0. 1-0.2/xHz pro- 
vided a good compromise, enabling the detection of closely 
spaced narrow peaks from long lived modes while keeping the 
number of spurious detections low around broad peaks arising 
from single short lived modes (Figure [Tlb-d)). 

We associated the frequencies, v, with radial or quadrupole 
(^ and 2) modes if e - 0.3 < i^/Ai/ mod 1 < e + 0.1, 
where e was derived using the /S.v-€ relation fitted to red giant 
branch stars by Cor s aro et a l, (2012) based on the formula- 
tion bv lMosser et al.l ( 1201 lbir (Figure[TTb-d). shaded regions). 

' http ://ai'chive. stsci . edu/kepler/ 
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Fig. 1 . — (a) Frequency spectram (noise nonnalized power) of a red giant 
branch star, with I'max and h.v indicated, (b) Small region of spectram of 
a red clump star. Spheiical degrees of the modes are indicated. Extracted 
frequencies (blue diamonds) are found as the peaks of the smoothed spectram 
(red curve). Black arrows show the location of the dipole modes apparent by 
visual inspection. All frequencies within the gray shaded regions are however 
assigned as radial and quadrapole modes in our automatic approach, (c) Same 
as panel (b), but for a red giant branch star of high inclination angle showing 
rotationally split dipole modes as close doublets, (d) Same as panel (c) but 
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The remaining frequencies were associated with dipole modes 
(I — 1), assuming negligible contamination from I = 3 modes 
(panel (e)). To ensure robust results, we discarded stars with 
fewer than five detected dipole modes overall. There were 
13 031 stars that passed this criterion. 

2.3. Measuring period spacings 

For each star, we calculated the pairwise period spacings, 
AP, between adjacent dipole modes, hence providing at least 
four spacings per star and usually many more. We plot these 
versus in Figure |2a). In this figure we show the AP 
regime that is populated by red giant branch (first ascent) stars 
(below about 100 s - see iBedding et an i201 1*). We find a 
remarkably clear upper envelope of the AP distribution in the 
range 6 < Ai///iHz < 16 at values of roughly AP= 70-90 s, 
which we interpret as a signature of the 'true' g-mode period 
spacing of dipole modes. 

To confirm this, we show the asymptotic g-mode period 
spacing for dipole modes, APg = ^/2-k'^{ J N /rdr)~^ , 
where N is the buoyancy frequency integrated over the radius, 
along a stellar model track representative for the typical stel- 
lar masses in our sample of stars (see Section |3). This is the 
period spacing of the g modes if they did not couple with the 
p modes. The alignment between APg from the model and 
the upper envelope of the observed AP is striking. The sharp 
edge of this upper envelope, and the fact that all tracks of 
representative masses fall a lmost on top of each other in this 
region (Section [I] see also iBedding et al.ll201 it IWhite et al] 
1201 U iMosser et al]|2012ah . is strong evidence that many of 
the detected dipole modes are essentially separated by the 
true g-mode spacing, which is now directly measurable for 
a significant fraction of stars. In order words, these modes 
resemble almost pure g-modes (see example in Figure [TJ)). 
Except for some special cases, we have previously only been 
able to infer the true g-mode spacing when making certain as- 
sumptions about the coupl ing between p- and g-modes (StelQ 
l20l2t lMosser"etani20 1 2ah . 

In addition, the over density of data points just below the 
model track (Figure |2^) gives us further confidence that we 
are essentially detecting the true g-mode spacing directly in a 
large number of stars. To illustrate this point, we show in Fig- 
ure[3]the period spacings of a sequence of / = 1 modes from a 
representative 1.5M0 model. The most common period spac- 
ings in the model are those close to the true g-mode period 
spacing, arising from the least bumped modes between the 
'dips' in AP where the slope of the AP{v) 'curve' is close to 
zero. The same effect could explain the slightly higher num- 
ber of modes spaced at about 50 s (as the slope of the AP{ v) 
curve also becomes zero at the bottom of each dip). These 
two extremes are marked as Max(AP) and Min(AP) in Fig- 
ure Etb). 

Stars at the bottom of the red giant branch with high val- 
ues of Av, towards the subgiants, typically show stronger 
coupling between p- and g-modes, and hence wider dips in 
the period s pacing, compared to what is shown in Figure |3] 
dPupret eTa l. 2009). In addition, their larger Av and AP 
makes their frequency spectra less dense, implying fewer 
modes spaced just below the true g-mode spacing. This seems 
to be supported by the absence of a clear over-density of 
points near the model track towards lower luminosity (higher 
Av) in Figure|2a). 

Another clear feature in this diagram is a large number of 
points with AP ^ 10 — 30 s. Many arise from rotationally 
spht modes and are discussed in Section [J!2l 



3. IDENTIFYING STELLAR POPULATIONS 

3.1. The AP-Av diagram 

To separate the stars into distinct groups, we need to as- 
sign a single representative value of AP to each star The 
median of all period spacings per star provided an efficient 
and robust measure for this purpose. The absolute values 
of AP from this method is most likely somewhat differ- 
ent than those found using previously publish ed methods 
(IBedding et al.ll201 lUMosser et al.ll201 laLl2012ah . but for our 
purpose that is not of concern. We show the result of the me- 
dian period spacing in Figure SJa), which for clarity only in- 
cludes modes that reach above six times the noise level (30 % 
of all stars). The masses indicated by color and shown as a 
histogram in the inset are derived from the scaling relation 

M/Mq ~ (!.,„ax/i^max.0)^(A!./Al^0)-4(T,ff/Teff,0)3/2, 

(2) 

which shows a relati ve scatter below 10% for samples 
of equal-mass stars dMiglio et al.l l2012h . This relation 
is derived from the scalin g relations for Av and z^max 
dKjeld^en & Bedding \\99^. For Teff, we used the values 
from|Pinsonneault et al. (2012). We found qualitative agree- 
ment for this plot when compa ring it with valu es derived using 
the method described in Moss eret al.l (1201 lal) . 

To compare qualitatively with our observations, we plot 
APg for a grid of stellar models in Fig ure Bib) using the 
MESA lM_pre_ms_to_wd test suite case dPaxton et al.l201 ll 
I2013h . In this dia gram the models evolve from right to left 
along the red giant branch. The hook at AP ~ 200 s along 
each track is the bottom of the red giant branch. When helium 
ignites at the tip of the red giant branch with Av values close 
to zero (outside the plotted range), the tracks quickly reap- 
pear on the helium-burning ZAMS (open diamonds; known 
as the zero-age horizontal branch for low-metallicity stars). 
Although the tracks were computed through this very rapid 
phase, for clarity we do not show this. About 2-3 stars per 
100 red clump stars a re expected to be in this transition phase 
(iBildsten et alJl2012h . Again, to simplify the plot the helium- 
core-buming phase is only shown from the helium-burning 
ZAMS until the core helium mass fraction has reached 0.14, 
except for the 1 .0 Mq track (dark blue), which shows the evo- 
lution all the way to the early phases of the asymptotic giant 
branch (dashed). As core convection due to helium burning 
ends, AP drops from 190 seconds to 100 seconds in less than 
1 Myr At the point where helium is completely exhausted in 
the centre of the star, which happens at AP ~ 80 s for the 
1.0 Mq track, the period spacing has been reduced back to 
something comparable to that of the red giant branch stars. 
We note that the change in AP as a function of mass along 
the helium-burning ZAMS follows quite tightly that of the 
helium core mass, which is lowest for M ^ 2.4 Mq without 
overshooting (IS. & M.II2013I) . 

From the speed at which models move across this diagram 
(Figurel^D), seen by the density of filled dots, we expect most 
observed stars to occupy the region, of 50 < AP< 90 (red 
giant branch stars with M < 1.8 Mq), 200 < AP< 250 (red 
clump stars with M < 1.8 M©), and 150 < AP< 200 (sec- 
ondary clump stars with M > 1.8 Mq). Although the median 
of our observed AP (panel a) is not exactly the same quan- 
tity as APg (panel b), the qualitative agreement between the 
models and the data is surprisingly good, and demonstrates 
that we can clearly separate the different stellar populations. 
How much the observed median AP deviate from APg de- 
pends on the fraction of individual period spacings that we 
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Fig. 2. — (a) AP versus Au for each pair of consecutive 1 = 1 modes (based on 0.1/jHz smootliing of tiie frequency spectra). Tlie black curve shows the 
asymptotic g-mode period spacing, APg, of a representative stellar model track below the red giant branch bump (see text), (b) Histogram of AP in the range 
8 < Ai///iHz < 12. 



O 
CD 

01 

a, 
< 



100 
80 

60 

40 

20 




: Model 


APg (1.5Ms) ■ 


n 




m 


n 



110 120 130 140 150 
Frequency (/^Hz) 



160 



Fig. 3. — AP versus frequency for each pair of consecutive 1 = 1 modes 
(diamonds) of a red giant branch model along the 1.5Mq track shown in 
Figure|2la) with Au = 10.8 /iHz and APg = 77 s. The dotted line marks 
the asymptotic period spacing. 



observe within, relative to outside, the dips (Figure |3), which 
again depends on both the strength of the coupling between 
the p- and g-modes and the number o f dipole modes per ra dial 
p-mode order that we detect (see e.g. lMosser et al.ll2012al) . 

One significant difference between panels (a) and (b) is the 
observed low-mass (blue and green) stars in the region brack- 
eted by AP 80 - 250 s and Ai^ ~ 5-7 /iHz, which is 
unexpected when compa red to the mode ls. These stars lie 
in the region where .Bildsten et al.l (1201 2 l) predicted the he- 
lium flashing stars to be while they establish stable core burn- 
ing. However, spot checks of the frequency spectra show that 
these stars are mainly red giant branch stars whose measured 
median AP are not representative of the actual spacing be- 
tween adjacent modes caused by our limited frequency res- 
olution. Hence, the median AP is artificially increased. To- 
wards lower values of Ai/, red giant branch stars show smaller 
period spacings, making it harder to resolve the individual 
mixed modes, which is somewhat exacerbated because we 
smooth the frequency spectra before we extract the frequen- 
cies. This also explains the lack of red giant branch stars with 

In Figure|4](c) we finally look at how metallicity affects the 
position in the AP-Ai^ diagram, in order to judge how a bias 
in metallicity could be introduced when selecting samples of 
stars from this diagram. It is quite clear that within the range 



of Z shown (corresponding roughly to —0.2 < [Fe/H] < 
0.2), the models do not show any significant shift in APg, 
and only a slight shift in Aiy, due to a change in mean stellar 
density. 

3.2. Internal rotation 

In Figure|5]we show the distribution of the observed median 
AP. Like Figure |2] it clearly shows a large fraction of stars 
with period spacings well below APg for red giant branch 
stars. Visual inspection of the frequency spectra shows that 
these are red giant branch stars with rotationally split dipole 
modes, as reported in a few st ars bv | Beck et akl (12 012) (see 
also the ensemble analysis by iMosser et al.l (l201 2b)). Fig- 
ures [Ttc) and (d) show some typical examples of such stars. 
Most of the pairwise period spacings for these stars are due 
to the rotational splitting of the dipole modes, which is why 
they show up with low values of the median AP. As illus- 
trated in Figure [lie) and (d), the dipole modes split into two 
or three peaks depending on their angle of inclination. We 
see two peaks when the inclination is high (close to 90 de- 
grees), three peaks are seen when the inclination is interme- 
diate (^40-60 degrees), and just one peak is se en when the 
star is viewed pole on (' Gizon & Solankill2003h . This sam- 
ple of rotating red giant branch stars provides an interesting 
selection for further studies into transport of angular momen- 
tum (Mosser et al. 2012b), age-rotation relations and the con- 
nection between surface and internal rotation including their 
angles of inclination. 

4. CONCLUSION 

We used a blind and fully automated approach to extract os- 
cillation frequencies for a very large sample (^ 13 000) of red 
giants observed by Kepler. From the extracted frequencies we 
measured period spacings between consecutive mixed modes 
of angular degree ^ = 1 . 

For a significant number of stars we were able to directly 
measure period spacings of the least bumped modes, suggest- 
ing that essentially pure g modes are now observable at the 
stellar surface. Because the mode Ufetime of the almost pure 
g modes is extremely long dBeck et al.ll2012l) . we can expect 
these modes will show up for an increasing number of stars as 
the time span of the Kepler observation increase. This indi- 
cates great prospects to further investigate the intricate details 
of stellar cores in red giants. 
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The median of the pairwise period spacings in each star 
from this simple but fast approach enabled an efficient way to 
distinguish red giant branch and red clump stars for the major- 
ity of stars in our sample. In addition, it enabled us to separate 
out a large fraction of stars showing rotational splittings. 

Our results demonstrate the great potential that this stel- 
lar sample holds for investigations into stellar structure and 
evolution, stellar populations of the Galaxy, and fundamental 
relations hips between stellar parameters such as age, r otation, 
0.10 K ' ' 
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Fig. 5. — Median period spacing of stars in FigureUa). Different popula- 
tions are indicated. The majority of stars in the range 80 < AP/s < 140 are 
red giant branch stars with 'incorrectly' detennined values of AP (see text). 
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Fig. 4. — (a) Median AP for each star versus its Au (based on 0.2/jHz 
smoothing of the frequency spectra). Red giant branch (RGB), red clump 
(RC), and secondaiy clump (2"'*RC) stars are indicated, (b) Curves show 
APg from theoretical models using MESA. Filled dots are separated by 
lOMyrs. Arrows indicate evolutionary direction. Masses of each track are in- 
dicated and their colors follow that of panel (a). The helium-burning ZAMS 
is marked by diamonds for each mass. The evolution from core to shell he- 
lium burning for the 1 MQmodel is shown by the blue dashed line, (c) Same 
as panel (b), but for 1 Mq tracks of three different values of the heavy element 
abundance. 



and metallicity. 
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